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We revise the spin-dependent neutralino-nucleus elastic scattering comparing the formalisms and 
approximations found in literature for the momentum transfer dependent structure functions. We 
argue that one of the normalized structure functions of Divari, Kosmas, Vergados and Skouras is 
all that one needs to correctly take into account the detailed nuclear physics information provided 
by shell-model calculations. The factorization of the particle physics degrees of freedom from the 
nuclear physics momentum dependent structure functions implied by this formalism allows for a 
better understanding of the so-called model independent method for setting upper limits. We 
further discuss the possibility of experiments with spin-dependent sensitivity like COUPP to test or 
set limits on the proton spin-dependent cross section in the framework of the stau co-annihilation 
region of the constrained minimal supersymmetric standard model. For this model with Aq = 0, we 
provide a fitting formula by which it is possible to convert an upper limit on the spin-independent 
cross section as a function of the neutralino mass directly into an exclusion plot in the (m 1( / 2 , tan/3) 
plane. 

PACS numbers: 95.35.+d, 12.60.Jv 



I. INTRODUCTION 

The nature of non-baryonic dark matter that seems 
to constitute the largest part of the matter in the Uni- 
verse is still unknown. If dark matter is formed by non- 
relativistic weakly interacting massive particles (WIMP) 
distributed in the halo of the galaxy, they should scatter 
elastically with the nuclei in a terrestrial detector [1]. A 
characteristic signal of the WIMP interaction is the pres- 
ence of an annual modulation in the event rate correlated 
with the motion of the Earth [2]. 

Experimental evidence of this modulation has been re- 
ported in the last years by the DAMA collaboration [3], 
and recently, also by the CoGENT collaboration [4] . The 
interpretation of these signals favors a light WIMP with 
mass around 10 GeV and a large spin-independent (SI) 
WIMP-nucleon cross section of order of 10~ 4 pb [5]. 
Other experiments, CDMS [6], XENON100 [7] and SIM- 
PLE [8], that anyway are not sensitive to the annual 
modulation, have reported upper limits that challenge 
the values of the cross section and mass statistically fa- 
vored by DAMA and CoGENT. 

If on the experimental side the situation is at least 
controversial [9, 10], on the theoretical side it is not 
less ambiguous. In the popular scheme of the minimal 
supersymmetric standard model (MSSM) with R-parity 
conservation where the lightest neutralino is a natural 
WIMP candidate, it is possible to accommodate a light 
neutralino with a cross section able to explain DAMA 
and CoGENT results while not contradicting other phe- 
nomenological constraints [11-14]. 

In supersymmetric models with unification conditions 
like the constrained MSSM (CMSSM) light neutralinos 
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with such a large spin-independent cross section are ex- 
cluded by other experimental constraints such as the LEP 
bound on the chargino mass. On the other hand, global 
fits that take into account accelerator, flavor physics and 
dark matter constraints, single out best fit points of the 
parameter space with a heavy neutralino [15-17]. 

In this paper we thus consider a region of the CMSSM 
parameter space, the so-called stau co-annihilation re- 
gion (tor)- In particular, we are interested to find out if 
present and future experiments can constrain this region 
by the spin-dependent (SD) elastic scattering. 

In the case of WIMP like the lightest neutralino (or 
any candidate with the same structure of coupling with 
nucleons), setting constrains on the SD couplings is, con- 
fronted with the SI case, complicated by the fact that: (a) 
there are two elementary cross sections, WIMP-proton 
and WIMP-neutron, that in principle should be con- 
strained at the same time and in a way that does not 
depend on the neutralino "composition" (the SI proton 
and neutron cross sections are to a very good approxima- 
tion equal) ; (b) in the formula for the neutralino-nucleus 
cross section the particle physics degrees of freedom are 
not factorized from the momentum dependent spin struc- 
ture functions (SSF) , thus when setting upper limits one 
is forced to fix the neutralino "composition" by the ratio 
of the couplings. Actually, problem (b) is at the root of 
problem (a). 

A solution to the problem (a) has been proposed in 
Ref. [18]. Thereafter the method has become the stan- 
dard way to derive limits on the SD WIMP-nucleon 
cross sections and to combine them from different ex- 
periments [19, 20]. 

We have discussed problem (b) in a previous paper [21] 
where it is evidenced that the foreseen factorization is ac- 
tually achieved by simply normalizing the standard struc- 
ture functions to their value at zero momentum transfer. 
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Here we show that the solution of problem (b) indeed 
gives a better understanding of the solution to the prob- 
lem (a) proposed by Rcf. [18]. In particular we show that 
there is no need of the assumptions made in Re. [18] that 
were object of criticisms Refs. [22, 23]. The method is not 
limited to the zero momentum transfer cross section but 
actually can incorporate the full momentum dependent 
structure functions. This is done in Section III. 

In Section II, and in the Appendix, we discuss various 
aspects of the momentum transfer dependent formalism 
and argue that some unnecessary complications of the 
standard formalism are at the origin of the aforemen- 
tioned problems. 

In Section IV we then discuss to what extent the limits 
on the single WIMP-nucleon cross sections derived by 
actual experiments like COUPP and XENNON100 can 
constraint the tcr- 

In Section V we give a parametrization of the 
SI neutralino-nucleon cross section in the stau-co- 
annihilation region that allows to translate an experi- 
mental upper limit into a bound in the (m 1( / 2 , tan/3) 
plane. 

Summary and conclusions are given in Section VI. In 
the Appendix we provide a detailed derivation of the for- 
mulas discussed in Section II. 



II. SD FORMALISM REVISED 
A. Structure functions and "form factors" 

Direct detection experiments employing odd nuclei 
with non-zero ground state angular momentum J, aim 
to constrain, in the case of absence of a positive signal, 
the spin-spin interaction of dark matter particles with 
the nucleons. Detailed nuclear shell-model calculations 
of the spin matrix elements in the zero momentum trans- 
fer limit (ZMTL), i.e. point-like nucleus, and of the SSF 
that account for the momentum transfer dependence re- 
sponse, have been carried out for many nuclei employed 
in actual experiments, see [24] for reviews. 

The differential neutralino-nucleus cross section, as 
function of the recoil energy of the nucleus E R — q 2 /2mA 
being q the modulus of the momentum transfer, has the 
general form 



dE R 



m A 
2^W 



* x A m\E R ). 



(1) 



Here A = SI or A = SD, tua is the mass of the nucleus 
with mass number A, \ia the neutralino-nucleus reduced 
mass and v the relative velocity. <7^)(0) arc the ZMTL 
total cross sections, to be discussed below. The function 
& X (E R ) accounts for the structure of the nucleus and is 
normalized to one in the ZMTL, $ A (0) = 1. 

For A = SI, <S> SI (E R ) = F 2 (E R ) where F{E R ) is the 
nuclear form factor. In Eq. (1), therefore, the nuclear 
physics is separated from the particle physics. 



For A = SD, in the standard formalism introduced by 
Engcl in Ref. [25], (see [24, 26, 27] for reviews), we have 



* S e D (Er) = 



S(E R ) 



(2) 



with 



S(E R ) = a 2 S 00 (E R ) + a oai S ai (E R ) + a 2 S n (E R ). (3) 

i, j = 0, 1 are isospin indexes and do and ai the isoscalar 
and isovector WIMP-nucleon scattering amplitudes writ- 
ten in the isospin basis. The ZMTL of the functions 
Sij{E R ) is Sij(0) ^ 1, they are not normalized to one 
and the function Soi for some nuclei can be negative. 
Particle physics and nuclear physics are not separated. 

These unpleasant features of the standard formalism 
are avoided with the formalism of Divari, Kosmas, Ver- 
gados and Skouras [28] . In this framework we can write 



$r(E R ) = 



HEr) 
J-(O) 



(4) 



with 



T{E R ) = a 2 F m {E R ) + 2a a 1 F 01 (E R ) + a?Fn (£?«). (5) 

Note that in this case -FV/(0) = 1 by construction. In 
Ref. [21] (see also [28-30]) we have remarked that the 
functions Fij(E R ) are practically identical in the recoil 
energy interval of interest for experiments, not only for 
light nuclei but also for medium-heavy and heavy nuclei, 



Foo(E R ) — F i(E R ) ~ Fu(Er). 
Thanks to Eq. (6), Eq. (5) reduces to 
^ D (E R ) = F n (E R ). 



(6) 



(7) 



Hence the SD "form factor" is determined by only one 
SSF. It does not depend anymore on the neutralino prop- 
erties as it happens in the SI scattering. 

The two formalisms are equivalent and connected by 



p (p , Sjj(E R ) 



(8) 



If the Sij are known also the Fij 



are known and vice 
versa. Eq. (7), anyway, allows for a drastic simplification 
of the formulas while retaining the exact informations of 
nuclear shell-model calculations. In literature, in spite of 
this, the formalism is largely overlooked. In some cases 
phcnomenological parametrizations are used. 

One example is the parametrization given in [31, 32] 



Fhs(qr n ) 



sin(qr 
qr 7 

047 



qr n < 2.55, qr n > 4.5, 
2.55 < qr n < 4.5, 



(9) 



with the nuclear radius r n ~ 1.0A 1 / 3 fm. 

Another example is furnished by the parametrization 
implemented in the code micrOMEGAs [33] for the case 
of nuclei for which the Sa are not available: 
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Figure 1. In blue line the normalized structure function Fu = Sn(q)/ Sn(0), for the nuclei 19 F, 73 Ge and I. The dashed 
line refers to the parametrization of Eq. (9) and the dashed-dotted line to the parametrization of Eq. (10). The variable in the 
abscissas is y = (qb/2) 2 , being q the momentum transfer and b = lfmyl 1 / 6 the oscillator size parameter. In the red abscissas 
the corresponding values of the recoil energies in keV are given. 



Fm0 = ^M = exi V^J' (10) 

where R A = 1.7A 1 / 3 - 0.28 - 0.78(yl 1 / 3 - 3.8 + [(A 1/3 - 
3.8) 2 +0.2] 1 / 2 ). These expressions are used also in recent 
literature [8, 34, 35] even in the case that the functions 
Sij or Fij are known. It is thus interesting to compare 
them with Fu. 

Figure 1 shows the normalized SSF Fu, Fls and F m o 
for one light nucleus, 19 F, one medium-heavy, 73 Ge, and 
one heavy nucleus, I, all of them largely employed in 
current experiments. The function Fu for 19 F is taken 
from Ref. [28], for 73 Ge is obtained from the function Su 
of Ref. [36], for 127 I from the function Su of Ref. [37] (set 
calculated with the Bonn A potential). 

In the abscissas we use the dimensionless variable 
V = (Qb/2) 2 where b = limA 1 / 6 is the oscillator size 
parameter. This variable is the natural one employed in 
shell-model calculations using harmonic oscillator wave 
functions. The functional form of Sij and Fij is typically 
a polynomial or a polynomial times an exponential in y 
or u = 2y [24, 28]. The recoil energy is easily found to be 
related to y by Er — 80 x y x A~ 4 / 3 MeV. For clearness 
we report also the corresponding recoil energies for each 
nucleus on a second abscissa. The interval < y < 1 cov- 
ers the recoil energies interval accessible experimentally 
but in the case of fluorine the relevant region is only up 
to y ~ 0.1. 



The approximation furnished by Fls is reasonable both 
at low recoil energies and at higher energies in the re- 
gion of the plateau, especially for the heavy nucleus. 
This is not surprising, for this parametrization was in- 
troduced [31] to fit the SSF in Xe and Nb [25, 38]. The 
approximation furnished by F m o is much worse in all the 
cases. A different Gaussian parametrization is given for 
example in Refs. [39, 40]. 

We stress again that for the nuclei for which the func- 
tions Fij or Sij have been published, there is no need 
of phenomenological fits or parametrization. The nor- 
malized function Fu accounts for the results of the most 
accurate spin structure function calculations and at the 
same time allows to separate the nuclear physics from 
the particle physics in SD the cross section. 

B. Differential and total event rate 

In SD scattering, given the neutralino-proton and 
neutralino-neutron cross sections = (/i 2 /7r)3|a Pjn | 2 , 
the total cross section at q = reads 

(ii) 

H p is WIMP-proton reduced mass and 



n p , n (0) =2\— -(S p ,„>, (12) 
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are the spin matrix elements of the proton and neu- parametrization proposed in Ref. [32] : 



tron groups. We remind that (S P) „) = (J, Mj 
J\S*JJ,Mj = J). In general both the SD WIMP- 
nucleon scattering amplitudes a p and a n (a Pi „ = (a ± 
ai)/2) and the nuclear matrix elements can have opposite 
sign, hence g = ±1 is the relative sign between |fi p (0)a p | 
and |fi„(0)a„|. An ab initio derivation of the SD cross 
sections using the formalism of Ref. [28] is given in the 
Appendix. 

In the SI case, for the neutralino we have cr^ 1 ~ <7„ = 
cr SI , the standard total cross section at q = is 



^(0) 



(AA 
Hp 



2 SI 



(13) 



The differential recoil rate is obtained by folding 
Eq. (1) with the velocity distribution function. We use 
the standard truncated Maxwellian [32]: 



v 



v v E 



f(v) 



(14) 



/(») = - 



cxp - 



(v - V E f 



exp 



(v + v E f 



k = y/irer{(z) — 2zexp(— z 2 ), z=-^-. 



^esc IS the escape velocity, v the velocity of the Sun v E 
the velocity of the Earth. 

Taking po — 0.3 GeV/cm 3 as the local dark matter 
density density, eo = ^^AVol^A/mA) the typical recoil 
energy and $ SI = F 2 (E R ), $ SD = Fn(E R ), we can write 



dR* 



-^(0) 



dt x 



dE R m x mA dE R 
dt x _ <f> x (E R ) 

»min(Bn) 



f> x (E R ) dv 
= / — /(«)• 

£ Jv^nfEn) V E 



dE R 

The total rate is simply given by 



m x rriA 
' E2 dt x 



J Ex 



(15) 
(16) 

(17) 
(18) 



F 2 (q) 



exp(- g 2 s 2 ), 



(19) 



. sinx COSX / 2 , 7 2 2 K 2 ( 

ji{x) = — 5 , r„ = \ c z + -TT 2 a 2 - 5s 2 fm, 

x z x V 3 

s = 0.9fm, a = 0.52fm, c = (1.23A 1/3 - 0.6) fm. 

In literature one can find other parametrization [41], or 
form factors obtained directly by shell-model calcula- 
tions [28, 29, 42]. We use here Eq. (19) because employed 
practically by all the experimental groups. 



III. MODEL INDEPENDENT UPPER LIMITS 

As an application of the previous formalism we dis- 
cuss the so-called model independent method for setting 
upper limits on neutralino cross sections and give an al- 
ternative proof of Eq. (13) of Ref. [18]. 

Let us consider a nucleus such that the SI rate is neg- 
ligible compared to SD one: in supersymmetric models 
the SD rate roughly dominate in nuclei with mass num- 
ber A < 20 while SI dominate at larger mass numbers 
due to A 2 proportionality [27]. We return on this point 
in the next section. 

We introduce the factors 



PqVq 
m x mA ' 

p A ^p,n(0) 



(20) 
(21) 



Eq. (17), with the aid of Eqs. (11), (20), (21), thus be- 
comes 



R 



SD 




(22) 



If an experiment with exposure £a = Ma x T, (Ma is 
the mass fraction of the element with mass number A 
and T the time of live data taking) have no statistically 
significant evidence, then an upper limit at some confi- 
dence level is put on the number of events N UL . For each 
unknown m x this is converted in an upper limit on the 
cross section requiring R x £ < N UL , that is 



The integration limits are v m { n (E R ) — vq i/E R /eo, 
"max = v esc , Ei = E th , E 2 = mm(E^ xp , E max ) where 
the maximal recoil energy is E maK = eo (v max /vo) 2 ■ The 
energy threshold E t h and E^ xp give the energy interval 
chosen by an experiment to analyze the data. For com- 
parison with a given experiment using specified nuclei 
and detection methods, if necessary, one should account 
in the previous formulas for the energy resolution and 
efficiencies that may depend on the energy. 

In the following we use the Helm form factor in the 




(23) 



The right-hand side of (23) is by definition the exper- 
imental upper limit on the neutralino-nucleus SD cross 
section, let us call <r A m as in Ref. [18], 



lim 
® A 



N 



UL 



b A t A D £A 



(24) 
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Furthermore, utilizing the same name of Ref. [18], we 
define the quantities 



lim(A) 
' p,n 



_lim 



{cry 



(25) 



Dividing both members of (23) by (24) and using the 
quantities (25) we arrive at 



< 1, 



(26) 



that is exactly Eq. (13) of Ref. [18] in the case of the 
allowed region in (<J p ,a n ) plane. 

In Ref. [18] the nucleon cross section limits in Eq. (25) 
are defined as basic quantities that then are combined to 
give Eq. (26). To do this it is necessary to assume that 
for a given nucleus it is possible to set separately limits 
on the SD-proton and SD-neutron cross sections even in 
the case that one contribution is clearly sub-dominant. 
These assumptions and the method were criticized in in 
Refs. [22-24]. 

In reality our derivation shows that such hypothesis 
are unnecessary and that the full justification of Eq. (26) 
only relies on the factorization of the particle physics 
from nuclear physics degrees of freedom and has a general 
validity 1 . 

Another common misunderstanding about Eq. (26) is 
that it is based on ZMTL total cross section and that it 
does not take into account the exact momentum depen- 
dent structure function. 

Actually, we see that using Fn, the correct behavior 
of the SSF can be taken into account in the upper limit 
CT iim(A) by the factQr f SD ; see Eq _ (24). 

On the other hand, the "upper limits" on the single 
proton or neutron cross sections, Eq. (25), are just useful 
quantities introduced to write Eq. (23) in the compact 
form (26). They become the actual experimental upper 
limits if, for the nucleus from which these are determined 
and in a specific WIMP model, one can prove that the 
protons contribution is dominant over the neutrons con- 
tribution or vice-versa (given the dominance of the SD 
rate over the SI rate). In general the exclusion curves on 
the single cross sections are fundamentally indicative of 
the experiments sensitivity and cannot constrain particle 
physics models in a universal way. 



IV. SD SCATTERING AND THE t C r 

To further clarify the last point, we choose a specific 
particle physics model, that is the constrained minimal 
supcrsymmetric standard model (CMSSM) with R-parity 



conservation. We consider the parameter space with 
fixed trilinear scalar coupling A$ — 0, positive Higgs mix- 
ing term (/x > 0) which is the benchmark supersymmet- 
ric theory for phenomenological and experimental stud- 
ies [43]. If the neutralino is required to furnish the cos- 
mological relic density inferred by WMAP [44] , then, for 
fixed tan/3 only specific regions in the (mi/2,mo) plane 
are left. In the (tcr) the lightest stau is almost degener- 
ate in mass with the neutralino and the co-annihilation 
of the two particles in the early Universe brings the value 
of the relic density in the favored WMAP interval. This 
parameter space is still untouched by direct detection 
experiments and LHC just started to explore it [15-17]; 
moreover it will be hard to probe it with indirect detec- 
tion methods such as 7-ray from neutralino annihilation 
in the halos [45-47]. 

The strips in the the plane (mi/2, mo) [43] for varying 
tan j3 from 10 to 50 in step of 5, are shown in the insert 
of Fig. 2(a). The strips and the cross sections are ob- 
tained with DarkSUSY [48], imposing WMAP constraints 
on the relic density 0.096 < flh 2 < 0.128, accelerator 
constraints on the lightest Higgs, m/, > 114 GeV and 
chargino mass m + > 103.5 GeV and the flavor physics 
constraint from bottom quark radiative transitions. 

In the same figure the SI neutralino-proton cross sec- 
tion as a function of the neutralino mass is shown. The 
SD neutralino-proton and neutralino-neutron cross sec- 
tion are shown in Fig. 2(b) and in the insert of Fig. 2(b), 
respectively. Two general features are worth noting: the 
SI cross section depends on tan /3 more strongly than the 
SD cross sections, the former varying by an order of mag- 
nitude and the latter by a factor less than 2; the SD are 
O(10 2 ) larger than the SI, in agreement with [57]. 

The neutralino field in the mass basis can be written as 
X° = N 11 B+N 12 W +N 13 H 1 +N 1 4H2, where N u are the 
elements of the matrix that diagonalizes the neutralino 
mass matrix, B, W° are the neutral gaugino fields and 
Hi, H 2 the neutral higgsino fields. In all the considered 
parameter space the neutralino is bino-like: we find nu- 
merically An ~ 0.99 > N 13 ~ 10~ 3 > Ai2,Ai 4 . This 
means that the coupling to the Z boson that is driven 
by the higgsinos couplings proportional to A13 and N14 
is heavily suppressed; the cross section is determined by 
squarks exchange. Analogously also in the SI case to 
the CP-even Higgs h and H are suppressed by A13 and 
A14 and the cross sections is thus mainly determined by 
squark exchange. Anyway, the couplings of the Higgs 
to down-type quark become (tan /3) 2 enhanced at large 
tan ft. The two contributions thus can be of the same 
order and SI cross section is more sensitive to variations 
of tan (3. 



3 F, 127 I and the t C r 



To discuss the relation between the SI and SD rates we 

1 This result was also implicitly obtained, with different notations consider the light nucleus 19 F that is known to furnish the 
and considering the case of general phases, in Ref. [30]. b est sensitivity to the proton SD CrOSS section [28, 49, 50] 
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Figure 2. (a) Spin-independent neutralino-nucleon cross section in the stau co-annihilation region with Ao — and tan /3 from 
10 to 50 in step of five as a function of the neutralino mass. The dashed line is obtained with the fitting formula, Eq. (29), 
with coefficients given in Table I. The insert shows for each value of tan/3 the strip in the plane (mu2, tUq) allowed by WMAP 
constraints on the relic density 0.096 < Qh 2 < 0.128 and satisfy accelerator constraints, (b) Spin-dependent neutralino-proton 
cross section and, in the insert, the spin-dependent neutralino-neutron cross section for tan j3 — 10, 30 and 50. 



and 127 I, that have both good SI and SD sensitivities. 

For 19 F we use the spin matrix elements of Ref. [28] 
that give ^ 9 (0) = 1.646 and f^ 9 (0) = -0.030. In 
this case the neutrons contribution in the SD rate can 
be safely neglected. We remark that the first nuclear 
shell-model calculation for 19 F [50] found (S p ) = 0.441, 
(S n ) = —0.109. The successive calculation of Ref. [28] us- 
ing a more realistic interaction, found (S p ) = 0.4751 and 
(S„) = —0.0087. The protons contribution is thus similar 
but the neutrons contribution is clearly negligible. The 
statement that the neutrons contribution is relevant, see 
for example [8] , in light of the more accurate calculation 
of Ref. [28], is doubtful. 

As reminded above, for light nuclei like fluorine the 
SD rate can be dominant over the SI, but this has to 
be checked in each particular WIMP model. We show 
the ratio i? SD /i? s/ for fluorine in Fig. 3(a). The SD 
rate is bigger by a factor up more than 2 at low and 
medium tan j3 but it is smaller than the SI rate at large 
tan /3; in any case the two rates are always of the same 
order of magnitude. The SI rate cannot be completely 
neglected at high tan/3 and for lower tan/3, neglecting it, 
one underestimates the total rate (see Ref. [51] for the 
case of general MSSM). The exclusion plots in the (m x , 
<7p D ) are inaccurate for the fca- In this case one has to 
draw an exclusion plot in the (<7p D , a SI ) plane for each 
fixed mass, the so-called mixed coupling approach [52]. 

Nuclear shell model calculations give for 127 I f2p 27 (0) = 
0.731 and r^ 27 (0) = 0.177 (spin matrix elements ob- 
tained with the potential Bonn- A from Ref. [37]). Al- 



though proton favoring, the neutrons group contribution 
to the nuclear spin is of the same order of magnitude. 
If the neutralino couplings to the proton and neutron 
are similar, the neutrons contribution to the nuclear spin 
must be considered. This indeed is what happens in the 
f C R where 0.75 < of^/of* < 0.9 [21] for tan/3 between 
10 and 50. Furthermore, a p < and a n > thus a 
cancellation in the SD rate is expected because the prod- 
ucts a p (S p ) and a n (S n ) are of the same order and have 
opposite sign. Fig. 3(c) shows the ratio of R ST> /R SI in 
127 I only considering the proton contribution, while in 
Fig. 3(b) both are included. Due to the A 2 proportional- 
ity, the SI rate always dominate by a factor from 4 to 25 
in Fig 3(c), but the cancellation makes the SD rate from 
2 to 3 orders of magnitude smaller than the SI, Fig 3(b). 

In the case of tqr, hence, iodine can only constrain the 
SI interaction. The exclusion plots in the planes (m x , 
(Jp B ), (m x , cr;p) or the combined (cp 50 , erf®) at fixed 
neutralino mass, derived using 127 I cannot constrain the 
^CRj for they are derived neglecting the dominant SI con- 
tribution or the equally important neutrons contribution 
that almost cancel the protons one. 



B. COUPP and the tcr 

The two nuclei discussed so far are the detecting 
medium of COUPP [53], a bubble chamber with CF3I. 
We use the latest data from Ref. [53]: an affective ex- 
posure of CF3I after cuts of £ — 28.1 kgxdays, 50% ef- 
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Figure 3. Ratio of the spin-dependent total event rate over the spin-independent rate R SU /R SI varying tan/3 in the stau 
co-annihilation region of the CMSSM with ^0 = and fi > 0. In panel (a) for 19 F. In panel (b), the ratio is plotted for 
127 Itaking into account both the protons and neutrons contribution in the spin-dependent rate; in panel (c) only the proton 
contribution is included. The points of the parameter space are the same as in Fig. 2. 



ficiency, E th = 21 keV, N UL = 6.7 at 90% confidence 
level and the same values of the velocities, wo = 230 
km/s, v esc — 650 km/s and an average velocity of the 
Earth ve = 244 km/s. In Fig. 4(a) the blue solid line 
is the present limit on SD WIMP-proton cross section 
derived from the fluorine fraction, while in Fig. 4(b) is 
the limit on the SI cross section derived from the iodine 
fraction. The blue dashed lines are limits extrapolated 
with the same N UL , 100% efficiency, effective exposure 
500 kgxyr and threshold at 7 keV. The red solid lines 
are the cross sections for tan f3 — 50, the orange ones for 
tan/3 = 10. 

The indication that we derive from Fig. 4 is that it will 
be unlikely for COUPP to probe the tcr by SD scatter- 
ing unless very large exposures of fluorine are achieved. 
On the other hand, a part of the parameter space will be 
probed by the SI scattering with iodine. This is not a lim- 
itation for the tcr since the two cross sections are clearly 
correlated. A constraint on <r SI automatically implies a 
constraint on er SD . As a matter of fact, should evidence 
be reported by more experiments and rates measured in 
fluorine, iodine and other elements like xenon, argon or 
germanium, the the full information on the SD sector can 
be reconstructed [21, 54]. 

As discussed in Section IV A, the SI rate in the 
tcr cannot be neglected for fluorine. Considering a 
mixed SI-SD approach with a fixed neutralino mass we 
obtain: 



SD lim(F) 



(c F ?t 



2+SD 



F a SI . 



(27) 



In analogy with Section II we set C = (/j,a/h p )A. At 
tan/3 = 50, where the SI rate is more important, the 
correction term on the right-hand side gets values larger 
than <7p D , in any case the largest values are of order 10~ 6 
pb. These values when compared with the present limits, 
crp im(F) ~ 10" 1 pb from COUPP and ~ 10" 2 pb form 
SIMPLE [8], are anyway negligible. Hence one should 
start to consider the SI rate only when the exposure is 
such that the sensitivity reaches the values of <7p D pre- 
dicted by the model. 



V. CONSTRAINING THE (m 1/2 , tan/3) PLANE 

In Fig. 4(b) we also show for comparison the present 
upper limit of XENON100 [7], which is the most stringent 
on the SI cross section. To be consistent with COUPP, we 
have calculated the XENON100 plot using the same val- 
ues of the velocities given above and the following data: 
effective exposure of 1471 kgxdays, energy threshold at 
8.4 keV and N UL = 5.62 at 90% confidence level deduced 
by the Feldman- Cousins method [56] with 3 events ob- 
served and mean background 1.8. As for the COUPP 
limits, we have calculated this curve using the total event 
rate without energy resolution function. Our curve dif- 
fers by few percent from published one, dotted line in 
Fig. 4(b). The latter is obtained with a statistical anal- 
ysis of the energy spectrum that take into account all 
the experimental uncertainties and with values of the ve- 
locities vq = 220 km/s, v esc — 544 km/s and ve = 232 
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Figure 4. (a) Spin-dependent neutralino-proton cross sections for tan/3 = 10 (orange line) and 50 (red line) in the stau co- 
annihilation region with Ao = 0, y, > 0. The blue solid line represent the present upper limit from COUPP, the blue-dashed 
line the extrapolated COUPP limit with an effective exposure of 500 kgx year and threshold at 7 keV. (b) The same as (a) 
but for the proton spin-independent cross section. The dashed-dotted line is the upper limit from XENON100 calculated as 
explained in the text, the dotted line the limit published in Ref. [7]. 



Table I. Coefficients for the fitting formula of Eq. (29). 



k 


(cr)fco (pb) 


(Pb) 


(<t)*2 (pb) 


(cr)fc 3 (pb) 


(a) k 4 (pb) 


2 


2.469xl0" a 


5.085xl0 -11 


5.432x10"^ 


1.783xl0" i3 


-6.089xl0~ ib 


3 


2.716xl0" 9 


-9.790xl0" 10 


3.92x10"" 


-6.413xl0~ 13 


-6.059xl0~ 15 


4 


1.395xl0~ 8 


-2.029xl0" 9 


2.143xl0" 10 


-6.711 xlO" 12 


9.481xl0" 13 



km/s. 

This exercise shows that for masses above 50 GeV the 
limits are more robust and less sensitive to the exper- 
imental details, statistical method to analyse the data 
and velocities (needless to say this is not true in the low 
mass region). In the high mass range m x > 50 GeV 
the exclusion limits are also robust against changes of 
the velocity distributions [55], being the major source of 
uncertainty a factor of two in pg. 

Since we have remarked above that the tqb. will be 
probably probed through the SI scattering, we further 
investigate what kind of information on the tcr parame- 
ter space can be extracted. We note from Fig. 2(a) that 
the SI cross section is a smooth decreasing function of 
the neutralino mass when mo and m 1 / 2 are varied along 
the WMAP allowed lines for fixed tan/3. Clearly it is 
also a continuous function of this parameter. Therefore 
we can look for a general fitting formula valid for all the 
values of tan/3. We first fit each a SI of Fig. 2(a) for a 
given value of tan /3 with the function 



A /100GeV\ fc 

2^* 



fc=2 



V 



(28) 



and than the coefficients (,k are fitted with a 4th order 
polynomial in tan/3. We thus find 



o- s/ (tan£,TO x ) = 



k=2 



A ( /100GeV\ fc ' 

.i=o \ m * / . 



(29) 

The coefficients of the fit are given in Table I. The 
fit obtained with Eq. (29) is shown in Fig. 2(a) with a 
dashed line. 

Analogously, the neutralino mass can be parametrized 
along the WMAP lines. We find that for all the values 
of tan /3 it holds 



0.44 mi/a - 15 GeV. 



(30) 



While the slope 0.44 is found for all the values, the con- 
stant negative term is an average value, since it presents 
a very mild dependence on tan /3 that anyway is not im- 
portant for what follows. Thus using Eq. (30) in Eq. (29) 
we end up with a formula cr s/ (tan/3, 7711/2) for the cross 
section in terms of the fundamental parameters mi/2 and 
tan j3. In last analysis, this allows to convert an upper 
limit on the event rate directly into an exclusion plot in 
the (mi/2, tan/3) plane. 
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Figure 5. Exclusion curves in the plane (mi/2,tan/3) for the 
CMSSM stau co-annihilation region with Ao — 0, fi > 0, set 
by an upper limit on the spin-independent neutralino nucleon 
cross section using the fitting formulas Eq. (29), (30). The 
regions to the left of the curves are excluded. The dotted line 
corresponds to the upper limit of COUPP(500 kgxyr), dashed 
blue line in Fig. (4)b. The other lines are extrapolations for 
XENON100 where the present effective exposure is multiplied 
by factors 5, 10, 100, 500. 



The result of this procedure is shown in Fig. 5, where 
the excluded regions are on the left of the curves. The 
COUPP upper limit with an effective exposure of 500 
kgxyear, dashed blue line in Fig. (4a), corresponds to 
the dotted black line in Fig. (5). The other curves are 
obtained for XENON100 considering an effective expo- 
sure to be 5, 10, 100, 500 times the present value of 1471 
kgxdays. The dot-dashed red line, corresponds roughly 
to the effective exposure of a future ton mass detector 
with 1 year operation and total acceptance cut of 40%. 
The extrapolation of COUPP, dotted black line, is ob- 
tained without any acceptance cut. 

We have to remark the limitations of the fitting for- 
mula. The coefficients in Table I have many particle 
physics uncertainties. First of all the cross section and 
the relic density were calculated with DarkSUSY with the 
default values of the hadronic matrix elements. Other 
codes can give slightly different values of the cross section 
for the same input parameters. Furthermore the depen- 
dence of SI cross section on not precisely known hadronic 
physics quantities can cause variations up to a factor five 
for a given point of the CMSSM parameter space [57]. 
There is a further dependence of the SI cross section on 
Aq. Anyway the choice Ao = is the benchmark case 
study also for direct searches of supersymmetric parti- 
cles at LHC: ATLAS and CMS typically present exclu- 
sion curves in the (mo, mi/2) plane with Ao — and fixed 
tan j3 [58, 59]. With the proposed formula, hence, one has 



a direct idea of the sensitivity of a direct detection ex- 
periment to one of the cosmologically favored region of 
the CMSSM parameter space in a complementary way to 
LHC. 



VI. SUMMARY 

In this paper we have reviewed the formalisms and the 
approximations found in literature for the treatment of 
the SD neutralino-nucleus elastic scattering. We argued 
that all that one needs to correctly take into account the 
detailed nuclear physics information provided by shell- 
model calculations is just one of the normalized structure 
functions of Rcf. [28] . 

We have shown that the factorization of the particle 
physics degrees of freedom from the nuclear physics mo- 
mentum dependent structure functions implied by this 
formalism allows for a straightforward proof of the gen- 
eral formula (26) proposed in Ref. [18] without the need 
of the assumptions that were criticized in Refs. [23]. 

We have further discussed the ability of some of the 
present experiments and their future upgrade to larger 
active masses (COUPP and XENON100) to constrain 
the stau co-annihilation region of the CMSSM. In this 
region of the parameter space the neutralino mass is in 
the interval 180-550 GeV and the SI cross section is a de- 
creasing function of the mass for 10 < tan/3 < 50, taking 
values in the range 10 -8 — 10~ 10 pb and it is still poorly 
constrained by experiments. The SD cross sections, with 
the proton and the neutron, are in the range 10~ 6 — 10~ 8 
pb. 

COUPP, although the high sensitivity of 19 F to the 
proton SD scattering and the fact that the SD neutralino- 
nucleon cross sections are larger than the SI neutralino- 
nucleon cross section, can constrain the model in its large 
mass phase only by the SI interaction with 127 I. The rea- 
sons are various: first, because of the A 2 scaling of the SI 
neutralino-nucleus cross section; second, there is a strong 
cancellation between the protons and neutrons contribu- 
tion in the SD neutralino- 127 I cross section; third, the 
active mass of 19 F is small. 

Furthermore, in 19 F, for the considered particle physics 
model, the SI rate is never negligible compared to the SD 
rate. In the case that the exposure were such that the 
model could be probed through SD scattering, an exclu- 
sion curve in the plane (m x , <7p D ) would be inaccurate. 

Finally we have given a fitting formula for the 
SI neutralino-nucleon cross section in the stau co- 
annihilation region as a function of the two fundamental 
parameters tan /3 and (10 < tan j3 < 50) that allows 
to directly convert an upper limit into an exclusion plot 
in the (tan/3, mi/2) plane for the case study Ao = 0. 
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Appendix: Proof of the formulas of Section II 

In this Appendix we derive the formulas given in Sec- 
tion II following the formalism of Ref. [28] but with a 
simplified and slightly different notation. 

The neutralino-nucleon SD cross section is determined 
by the axial part of effective lagrangian. At the nucleon 
level, in the isospin representation that is convenient for 
nuclear physics calculations, we can write 



(Al) 



The operator f 3 act as r 3 \p) — \p), f 3 |n) = — \n) and 
1 is the identity operator in isospin space. Thus for 
N = p,n the isospin operator gives (ao ± &i)/2 = 
a P ,n = J2 q d q Aq( N ^ where d q is the effective coupling 
with quarks and Aq( N *> the spin fractions of the nucleon 
carried by the quarks. We do not discuss further the 
physics involved at the nucleon level, see Refs. [49, 57]. 

Taking the non-relativistic limit we get the neutralino- 
nucleus spin-spin interaction 

A 1 

V = 4§ x ' Yl o ( flo1 + a i f ')Mr - r< ). (A.2) 



Here Si and are the spin and coordinates of the i-th 
nucleon. In literature sometimes factors Gf/V%, GfV% 
or are extracted from oo,i. To simplify the for- 

mulas, we adopt instead the convention that all the cou- 
plings are included in a .i. 

The spin operator of the neutralino operates on eigen- 
states of the spin \s), while conventionally all the an- 
gular momentum operators of the nucleus are evaluated 
in state with the maximal value of the z projection, 
(6) ee (J,Mj = J\O z \J,Mj = J). The nuclear wave 
function depends also on the isospin and the coordinates 
of the nuclcons 



\A) = \ J,M J = J,T 3 ,r 1 ...r A ). 



(A.3) 



The elastic differential cross section in the center of 
mass frame and the total cross section, in the case that 



there is no angular dependence of the amplitude, are 
given by 



where the scattering matrix element, with \A, x) — 
\A)\s), is 



M 



(Ax\ J 



dve-^ r V\A, x ). (A.5) 



For two spin operators acting on different spaces, the 
average over the initial directions of modulus squared of 
the scalar product is |S a • Sb| 2 = ^S^S 2 , hence 



1 



\M\^-Msl) s {S ) A , 



(A.6) 



We have defined the operator 

A 

2"' '' "~" 2 



(A.7) 



with 



n = ^ lS ie - jq r * , Oi - ^S je - Jq r * . (A.8) 



To evaluate (S ) a we note that for a vector operator, 
the matrix elements in states | J, Mj = J) are related to 
the reduced matrix elements by [60] 



(J,J|OV,J) = s7 ^ T IVI|0||J)| 1 
It follows: 



(A.9) 



(J, j|6 2 | j, j) = 1±1\(j,j\6 z \j,j)\ 2 . (Aio) 

We thus define the momentum dependent matrix ele- 
ments 



J 



From Eqs. (A.7)-(A.ll) we find 

(S 2 ) A = ^\a o n o (q) + a^iq)] 2 . 



(All) 



(A.12) 



Obviously, (s 2 ) s = s (s + 1) = 3/4. Eq. (A.6) thus takes 
the form 



|M| 2 = -|aof>o(<z)+aitti(<z)| 2 - 



(A.13) 



11 



Expanding the square and factoring out the zero mo- 
mentum values, we introduce the normalized structure 
functions Fij(q): 



Cli(q)Clj(q) 
fii(0)fij(0)' 



and find 



Furthermore, using ao,i = a p ± a„ and Eq. (12), 
Eq. (A. 16) becomes: 



\ M \2 = 4 ^_p (ap ( Sp ) + a „(S„)) 2 Fu(g). (A.19) 



(A.14) 

For a single nucleon (A.19) reduces to 3|a Pi „| 2 , hence 



SD 



3 — \a v 



(A.20) 



\M\ 2 - -(a 2 Q 2 (0)F 0Q (q) + 2a a 1 Q (0)Q 1 (0)F 01 (q) 



a 2 Q 2 (0)F n (g)). 



(A.15) 



By reason of Eq. (6), we can make the approximation 



\M\ 2 ~ -(a O (0) + oifii(0)) J Fn( g ). (A.16) 



Taking q = in Eqs. (A. 11) and using f 3 |p) = +|p) and 
f 3 |n) = — |n) we have 



<)„,„, ,^/i+l(£ i isf) A = 2 ] l^±l((S l ,) -«S„» 



= O p (0)+O„(0), 

l — l 

= n p (o)-n„(o). 



Eq. (12) is thus proved. 



(A.17) 



2 A /^±i«S p )-<S n » 



(A.18) 



Eq. (11) follows from Eqs. (A.19)-(A.20). 
Finally, with the substitution dfl = 



4tt 



■dq 2 



^^-dE R in Eq. (A.4) and using Eqs. (A.19) and (11) 

we obtain Eq. (1). 

The present formalism and the standard formalisms 
are equivalent and connected by 



Sij{q) 



2J+ 1 
(1 + <%)8tt 



(A.21) 



For a given nuclear wave function they furnish the same 
cross section. 

In last analysis, the difference resides in the way by 
which the multipole decomposition in vector spherical 
harmonics of the operator (A. 7) is carried out. In the 
standard formalism this done in terms of the operator 
T£ l 5 and C\, see Refs. [24-26] and references therein for 
explicit formulas and meaning. Both the operators con- 
tain the couplings ao and a\, thus the modulus squared of 
each contains terms proportional to a 2 ,, a 2 and the inter- 
ference aofti. The function S(q) in terms of the functions 
Sij(q) arises after a rearrangement of these terms. 

In the formalism of Rcf. [28] the multipole decompo- 
sition in vector spherical harmonics is done on the oper- 
ators (A. 8). Anyway, keeping separated the terms in a 
and on has the advantage that the squared of the ampli- 
tude is always a perfect square, see Eq. (A. 13), and limit 
q = is reached in a more transparent way. 
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